Theoretical prediction of rotating magnon wavepacket in ferromagnets 
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We theoretically show that the magnon wavepacket has a rotational motion in two ways; a self- 
rotation and a motion along the boundary of the sample (edge current). They are similar to cyclotron 
motion of electrons, but unlike electrons the magnons have no charge and the rotation is not due 
to Lorenz force. These rotational motions are caused by the Berry phase in momentum space from 
magnon band structure. Furthermore, these rotational motions of the magnon give an additional 
correction term to the magnon Hall effect. We also discuss the Berry curvature effect in the classical 
limit of long- wavelength magnetostatic spin waves having macroscopic coherence length. 

PACS numbers: 85.75.-d, 66.70.-f, 75.30.-m, 75.47.-m 



Introduction — Spin wave (magnon) in an insulating 
magnet is a low-energy collective excitation [l|, Q- It 
has recently been focused as a tool for the spintronics 
application, because it can have a good coherence, com- 
pared with the spin current in metals. The motions of 
the magnons are now measurable in a time- and space- 
resolved way with reasonable accuracy. It can be experi- 
mentally generated and detected via the spin Hall effect 
0, and for the spintronics application, a precise spatial 
and temporal control of spin wave is desired. 

In the present paper, we theoretically find that the 
motion of magnon wavepackets in insulating magnets un- 
dergoes a self-rotational motion and a rotational motion 
along the edge of the sample (Fig. [Ija)). The latter mo- 
tion gives rise to the thermal Hall effect of magnons ■ 
This phenomenon is due to the Berry curvature in mo- 
mentum space, representing the topological structure in 
the magnon bands. We theoretically predict that if a 
magnon wavepacket is excited in the vicinity of the edge 
of a sample, it will move along the edge. It is expected to 
be visible in some magnets having macroscopic coherence 
length (~ 10mm) Q, and it should be a powerful tool for 
exploring effects of Berry phase in momentum space. 

We calculate the transverse thermal transport coeffi- 
cients for a magnon system by two methods: the semi- 
classical theory and the linear response theory, by anal- 
ogy with an electron system. We show that both theories 
give the same result and the thermal Hall conductivity 
can be written in terms of the Berry curvature. Our 
theory includes the contribution of magnon rotational 
motion, which has been overlooked in the previous the- 
ory of magnon thermal Hall effect 4, sj . From this we 
show that the thermal Hall effect of the magnon arises 
from the edge current of the magnon. We apply our 
theory to various magnons, including both the exchange 
spin wave (quantum-mechanical magnon) e.g. in a fer- 
romagnet LU2V2O7, and the magnetostatic spin wave in 
yttrium-iron-garnet (YIG) films. The results for ther- 
mal Hall conductivity in LU2V2O7 roughly reproduces 
the experiment [5l|. Throughout this paper we consider 
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FIG. 1: (Color online) (a) Self-rotation of a magnon 
wavepacket and a magnon edge current, (b) The magnon 
near the boundary proceeds along the boundary, irrespective 
of the edge shape, (c) Magnon edge current in equilibrium, 
(d) Under the temperature gradient, the amount of the trans- 
verse heat current are not balanced between the two edges, 
and a finite thermal Hall current will appear. 



localized spin systems on a two-dimensional lattice, and 
assume that there is no interaction between magnons. 

Semiclassical theory — The dynamics of a wavepacket 
of electrons in a periodic system can be described by the 
semiclassical equation, including the topological Berry 
phase effect [8]. When a force is exerted onto the elec- 
tron, there occur various intrinsic Hall effects due to the 
Berry phase. In analogy with this, we construct the semi- 
classical equation of motion of magnons. We consider a 
wavepacket of a magnon which is localized both in real- 
and momentum-space. If there exists a slowly-varying 
potential U (r) for the magnons, they feel a force. Follow- 
ing I2I , we derive the semiclassical equations of motion 
for the magnon wavepacket as: 
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k X f2„(fc), hk ^ -VUir), (1) 



where fln{k) is the Berry curvature: 11„(A;) = i (^^| x 
1^^) with \un{k)) being the periodic part of the Bloch 
waves in the nth band and Enk is the nth band energy 
of magnon. For the potential U{r) we cannot use the 
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electric field because the magnons have no charge. In- 
stead we focus on the boundary of the system, which can 
be regarded as a confining potential U{r). On the edge 
along y-direction, for example, the gradient of the con- 
fining potential dxU{r) produces an anomalous velocity 
k X r2„ = —hr^ dxU {r)Q.n^z{k)y in Eq. ([1]). By summing 
it over all the occupied states, we get an edge current 
/, = -jldxdxU(r) {^EnMPi^nk + Uir))n^^,{k)) = 

~wJ2n,kI^^'^^Pi^)^n,zik), where x ^ a and x ^ b 
represent the inside and the outside of the system, re- 
spectively, such that U{a) = and U{b) = oo, p{e) is 
the Bose distribution function p{e) — (e^^"^)/''^^ — 1)"^, 
/cb is the Boltzmann constant, fi is the chemical poten- 
tial and T is the temperature. Since the magnon edge 
current is independent of the edge direction, we simply 
write this as /. This magnon edge current / is indepen- 
dent of the form of the confining potential. Therefore, the 
edge current circulates along the whole edge (Fig.[TIb)). 
Strictly speaking, this approach is applied only when U 
is slowly varying. Nevertheless, as / does not depend 
on the form of U{r), this remains valid even when U{r) 
is representing a hard wall and rapidly varying. Such 
kind of approach has been successful in quantum Hall 
systems tlOj. We can therefore expect that the similar 
approach is successful also for magnons; we can use the 
slowly varying confining potential which is much easier 
for theory in order to predict physical phenomena occur- 
ing irrespective of the details of the potential. We note 
that if the coherence length is short, the effective system 
size is given by the coherence length of the magnon, and 
the edge current is also confined within this length scale. 

If either /i or T varies spatially, this otherwise cir- 
culating current will no longer cancel in the interior of 
the system. This causes the thermal Hall effect as we 
show in the following. The magnon current and en- 
ergy current due to the edge current can then be writ- 
ten respectively as j = V x -f^Y.n,kI^^ p{£)^n{k)d£, 
jE = V X j^Y.n,k I^k'^Pi^)^ri{k)de. Various ther- 
mal coefficients are derived from these equations. If 
the chemical potential fi or the temperature T spatially 
varies, the magnon current and heat current is induced 
via the Bose distribution function. For instance, in the 
presence of a temperature gradient in the y direction, 
the magnon current in the x direction arc written as 

Other thermal transport coefficients can be obtained in 
the same way. Now we write the linear response of the 
magnon current and the heat current as 

j - Lu [-VU - V/i] + Li2 [TV {^)] , (2) 
Jq =Li2[-Vi7- Vm]+L22 [rv(^)] , (3) 

where the heat current jq is defined as jq s j'e — fij. 
We take p = here, because the magnon number is not 



conserved. The transport coefficients can be written as 
L-]'^-^-^^T.^nAk)cM, (4) 

where V is the area of the system, p„ = /9(e„(A;)), Cq{p) — 
{log{l + t-^)y dt,q = i+ j ~2 with i,j = 1,2. For ex- 
ample, co(/o) = p, ci(p) = (1 + p) log (I + p)-plogp, and 

C2(p) = (1 + p) (log - (logp)' - 2Li2(-p), where 

Li2(2;) is the polylogarithm function. Thus the thermal 
Hall conductivity = ^22 /"^ '^^^ be obtained as 

«"^ = -^E^2(Pn)^^„,.(fc). (5) 

From Eq. ([5|) we can see that comes from the Berry 
curvature n„{k) in momentum space. Therefore, if the 
energy band is close to each other, i.e. near the band 
crossing, this gives a large contribution to k^^. 

Thus, the thermal Hall conductivity solely comes from 
the edge magnon current. In equilibrium (Fig. [Ijc)), 
there exists the edge current of the magnon due to the 
confining potential gradient. This magnon edge current 
circulates along the boundary, giving no net thermal cur- 
rent across the magnet. If the temperature gradient is 
applied (Fig. [Hd)), the balance of the contributions to 
the heat current from the two opposite edges will be bro- 
ken, and finite thermal Hall current will appear. Because 
the temperature gradient is a statistical force, it neither 
exerts a force to the magnons nor deflects the wavepacket 
in the bulk. The edge current should be observed experi- 
mentally by using time and space resolved Brillouin light 
scattering technique pj. 

Linear response theory — There is a discrepancy be- 
tween Eq. (0) and the formula obtained in 0, @ using 
the Kubo formula. In 0, [Hj , the term from magnon rota- 
tional motion is missing, as we discuss in the following. 

Because the temperature gradient is not a dynamical 
force directly acting onto the particle, but a statistical 
force, its theoretical treatment requires some care. In 
the linear response theory for electronic systems includ- 
ing heat current and temperature gradient, it is conve- 
nient to introduce a fictitious gravitational field ip [Hj, 
which exerts a force to the wavepacket, proportional to 
its energy. Such formalism can be applied to the magnon 
system, taking into consideration several differences that 
the magnon has no charge and is not a fermion but a bo- 
son. As a result, the transport coefficients for magnons 
consist of two terms: one term from a deviation of a par- 
ticle density operator from an equilibrium calculated by 
the Kubo formula, and the other term arising from the 
deviations of current operators which are linear in ap- 
plied fields. The latter term is expressed in term of the 
reduced orbital angular momentum of magnons ^ (rxv), 
but was missing in 0, [Hj . This latter term has been dis- 
cussed in the context of orbital magnetization in electron 



3 



systems 12h17| . The calculation is in parallel with that 
for electrons |l2l[l3| . and its details will be presented else- 
where. As a result, the thermal Hall conductivity from 
the linear response theory is identical with that from the 
semiclassical theory in Eq. ([S]). 

The new term to the linear response theory corre- 
sponds to the reduced orbital angular momentum of 
magnon. It consists of two parts: the edge current and 
the self-rotation of the wavepacket. The reduced angular 
momentum for the edge current per unit volume is 
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and that for the self-rotation is calculated in analogy with 
the electron system [3| as, 



{H — £nk) 



dUn 



(7) 



Namely, in addition to this edge current, we find that 
the magnon wavepacket rotates around itself and induces 
orbital angular momentum, because of the Berry phase in 
momentum space. Thus the magnon in equilibrium has 
in general a nonvanishing orbital angular momentum due 
to the Berry curvature. This magnon orbital motion can 
be regarded as a generalized cyclotron motion, whereas 
the magnon feels no Lorentz force and cannot have a 
cyclotron motion in the same sense as that of electrons. 
In this respect, this motion is purely due to the magnon 
band structure. This effect is common in various wave 
phenomena like electrons 9], photons [l^, and so forth. 

Orbital motions of electrons give rise to a magnetic 
moment due to the electron charge. On the other hand, 
magnons have no charge, but have a magnetic moment. 
Because the magnon carries magnetic dipole, the rotating 
magnon wavepacket can be regarded as a circulating spin 
current. Hence as is similar to the spin Hall effect, and 
its insulator counterpart, i.e. the magnetoelectric effect 
in noncollinear spin structure (2^, the rotating magnon 
wavepacket should accompany a polarization charge. It 
requires the spin-orbit interaction, i.e., Dzyaloshinskii- 
Moriya (DM) interaction. This effect is dual to the rota- 
tion of electric charge, producing a magnetic dipole. 

Thermal Hall effect of Lu2 V2 O7 — Now we apply 
our results to the ferromagnetic Mott- insulator LU2V2O7 
with a pyrochlore structure. This material has spin-1/2 
V'*''' ions with the DM interaction. The collinear ferro- 
magnetic ground state is stable because the total DM 
vectors of the bond sharing the same site is zero [Sj], 
and the effective spin-wave Hamiltonian is written as 
Hcs = -JSi- Sj+Dij ■ {Si X Sj)-g^iBH ■Y.i St, 

where denotes the nearest neighbor pairs, J is the 
exchange interaction, D is the DM vector, g is the g- 
factor, /xb is Bohr magneton, and H is the magnetic 
field in the z direction. We focus on the tempera- 
ture regime much lower than the Curie temperature 
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FIG. 2: (Color online) Numerical results, (a) thermal Hall 
conductivity of the magnon in LU2V2O7 in the magnetic field 
at T = 20[K]. (b) Berry curvature n^{k) with n = 0, 1, 2 for 
the magnetostatic spin waves in YIG. Eigenmode dispersions 
are shown in the inset. 



Tc — 70 [K], for existence of well-defined Bloch waves 
of magnons. We can then assume that the contribution 
from the lowest band is dominant, whose Berry curva- 
ture is rii.z ~ ~'^/2T^^^^ + + 2^2) as calculated in 
5], with A being a quarter of the lattice constant. Us- 
ing this, we can estimate the orbital angular momentum 
of the magnon from both the self-rotation motion L|°'^ 
and the edge current L°'^s°. Near k = 0, the lowest-band 
dispersion is quadratic and we can introduce the effec- 
tive mass of the magnon of the lowest band ml , defined 
as m* = lT?{d^enk/dk^)^^. The orbital angular momen- 
tum from the self-rotation motion of the magnon per unit 
volume is calculated analytically. 



JSm\D 1 (k^T 



5/2 



Lis 



hA J 327r3/2 \^ js 

' (8) 

We obtain " ~ -0.009ft and Lf^" ~ -hO.OOSftper unit 
cell. We can also calculate the thermal Hall conductivity 
, assuming that the contribution of the lowest band is 
dominant. The resulting is shown in Fig.l^Ta), which 
roughly agrees with the experimental results [5|. 

Magnetostatic spin waves — Our theory is also ap- 
plied to the magnetostatic spin waves in a ferromagnet. 
In the magnetostatic spin waves, the wavelength are suf- 
ficiently long and magnetic anisotropy is mostly deter- 
mined by the demagnetizing field which is dependent on 
the sample shape. The exchange coupling is negligible 
because of the long wavelength. This anisotropy due 
to the demagnetizing field plays the similar role as the 
spin-orbit coupling, and gives rise to the Berry curva- 
ture contributing to the Hall effect of spin waves. As 
an example, we consider YIG films, magnetized to the 
saturation in an arbitrary direction by an external mag- 
netic field. YIG is a ferrimagnetic insulator, and the 
spin wave in YIG can propagate over centimetres. The 
spin wave mode is expressed as a plane wave: m.(r,t) — 
m(z)exp(i(fc • ry — ujt)), where m is a two-dimensional 
vector perpendicular to the saturation magnetization 
Afo, 2: is a coordinate perpendicular to the film, ru is 
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the coordinate within the film, and a; is a frequency of 
the spin wave. The linearized Landau-Lifshitz equation, 
coupled with the Maxwell equation with boundary con- 
ditions for the film, is cast into the integral equation 2l| : 

ujH'm{z) — ojm /_f'^y2 ^■^'^(-^' -^')"^(-^') = (^o-ym{z). Here, 
u!h = jHf), lom — jMq, L is the film thickness, G{z,z') 
is 2 X 2 complex matrix of the Green's function defined in 
2l| . Uy is the Pauli matrix, 7 is the gyromagnetic ratio. 



Hq is a static magnetic field. 

This integral equation is a generalized eigenvalue prob- 
lem due to the presence of ay. The calculation of the 
Berry curvature then requires some modifications. Fol- 
lowing the prescription of the wavepacket dynamics , 
we introduce the Berry curvature of the magnetostatic 
spin wave Vtlik) = -ea/3^Im (^^^^rf- cFy ^5^) ' where 
ea/3-y is the totally antisymmetric tensor, n is a band in- 
dex of the spin wave, and the bra-ket product refers to 
a usual inner product of vectors and an integral over z. 
In some cases, this Berry curvature vanishes because of 
symmetry. This occurs when Mq lies in the plane. In this 
case the system is invariant under product of the time- 
reversal operation and the n rotation within the film, and 
this symmetry forces the Berry curvature rin(fc) to be 
zero. On the other hand, when Mo is not in the plane, 
the Berry curvature is expected to be nonzero for any 
modes, leading to the rotational motion of the wavepack- 
ets and the Hall effect. Actually, if Mq is perpendicular 
to the plane, i.e., for the magnetostatic forward volume 
wave (MSFVW) mode, we can calculate the Berry curva- 
ture. From Ref. p3| . the frequency w = a;„ for nth eigen- 
mode is determined by y/pta,n{ '^^^ _(_ i^) = 1 where 

_ i^Af^g — I Qj^^^ n = 0, 1, 2, • • • , and is shown in the 
inset of Fig. [DJb) for n 0, 1, 2. The modes are 

m„fe(z) = V - X fe))i (9) 



where TV is the normalization constant determined by 
(''T^n.fcl <yy \'mn,k) = 1- Wc then obtain for n-th mode; 
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fl^ is evaluated numerically, and the results are shown in 
Fig.HJb). We have thus confirmed that the Berry curva- 
ture is nonzero for the MSFVW mode, and the rotations 
of magnon wavepacket are predicted to occur. 

Conclusions — To summarize, we found that the 
magnon wavepacket undergoes rotational motions in two 
ways: self-rotation and motion along the edge. The 
latter is responsible for the magnon thermal Hall ef- 
fect. These rotational motions are due to the Berry 
phase, and is similar to the electron cyclotron motion, but 
without a Lorentz force. The present theory is applied 
both to the exchange spin wave (quantum-mechanical 



magnon) e.g. in LU2V2O7, and to the classical magne- 
tostatic waves e.g. in YIG. This effect is expected to be 
observed via space- and time-resolved observation of the 
magnon wavepacket with macroscopic coherence length. 

Similarly to anomalous Hall or spin Hall effect, the 
Berry curvature in the magnon systems is enhanced near 
band crossings, where the magnon frequency in a focused 
band is close to that of other bands. One can design the 
magnonic crystals jl^l so that the magnon bands have 
a band crossing. One can then expect to have a promi- 
nent rotational motions of magnon wavepackets, if its 
wavenumber is close to the band crossing. 
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